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Целью работы являлось создание метода оценки и прогнози- 
рования надежности сложных систем на основе статистиче- 
ского имитационного моделирования их работы. Использован 
метод рисков с учетом потерянной выгоды. Предложены зави- 
симости для определения потерь от отказов при функциони- 
ровании системы в зависимости от степени ее безотказности. 
Разработана программа для проведения имитационного моде- 
лирования с учетом рисков от надежности — "Вена Шу". 
Проверка разработанного метода проводилась на основе ими- 
тационного моделирования работы измельчителя компонен- 
тов комбикормов типа У17-УКИ. Полученные результаты 
отражают реальные процессы, происходящие при эксплуата- 
ции технических объектов. Метод полностью формализован и 
может быть использован для оценки надежности любых тех- 
нических систем. 


Ключевые слова: надежность, риск, статистическое имита- 
ционное моделирование. 
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Введение. Отказ технической системы неизбежно ведет к потерям: производство останавливается или сокра- 
щается, отказавшая система требует ремонта, а последствия аварий —Щ_ ликвидации. Кроме того, эксплуатация нена- 
дежной техники может оказывать негативное влияние на окружающую среду и людей. 

Экономическая оценка от ненадежности оборудования в настоящее время проводится в два этапа. Сначала 
рассчитываются показатели надежности системы или машины, а затем определяются убытки. Такая процедура до- 
вольно длительная и трудоемкая. Если использовать такое понятие как риск, то экономическую оценку ненадежности 
можно проводить быстрее и точнее. Риск является неизбежным атрибутом эксплуатации техники. Он является одним 
из важнейших показателей безопасности. Риск, возникающий в результате отказов техники, называется техногенным. 

Риск можно рассматривать как вероятность некоторых неблагоприятных событий, но обычно под риском по- 
нимают оценку ожидаемого вреда (потерь) от неблагоприятных событий. В большинстве случаев риск оценивается 
денежными единицами, хотя могут быть и иные случаи. 


Степень риска может быть оценена следующим образом [1-3] 
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в=0(1)-С=[1-Р(Е)]С 


В(Е)=Уус,-РИЕ) (1) 
где ^ — величина риска, Р(!) — вероятность безотказной работы системы, О(1) — вероятность появления отказа, С — 
величина потерь, возникающих в результате отказа, с; — потери от отказа ]-ого элемента, р;#) — вероятность безот- 
казной работы /-ого элемента. 

Постановка задачи. Вопрос заключается в определении О(1) или Р(1). Чаще всего определяется Р(!). Для это- 
го строится структурная иерархическая схема системы (машины), состоящая из конечного количества элементов и на 
основе Р.(#) (вероятность безотказной работы одного элемента) определяется Р(!). 

Основные формулы для различных соединений. Элементы в системе могут быть соединены последова- 
тельно, параллельно или смешано. Но, в конечном итоге, это будет сложное последовательное соединение. При таком 
соединении вероятность безотказной работы системы равна произведению вероятностей безотказной работы всех 
элементов [4, 5]. 


Ре(1) = Пар (Е) (2) 


где р:()— вероятность безотказной работы. 
При параллельном соединении т элементов вероятность безотказной работы: 


ре(#)=1-П;1-РЕ) (3) 


Поэтому параллельное соединение используется для эффективного резервирования системы. Для смешанного 
соединения (п последовательных групп из т параллельных элементов) вероятность безотказной работы: 


== п п 
ре(г)= Пи Пи -рх#}} (4) 
Каждый элемент характеризуется тремя показателями: наработкой на отказ Т,, временем восстановления Г,и 
ценой. Т, и Т, случайные величины, которые характеризуются соответствующими распределениями. 
Программная реализация процесса. Решение задачи в такой постановке лучше всего проводить методом 
статистического моделирования [4,5]. 


Для оценки надежности системы и учета рисков была разработана программа для имитационного моделиро- 
вания ее работы. Блок-схема программы приведена на рисунке 1 [6]. 
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Рис. 1. Блок-схема программного комплекса КепабИиу 
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Достоинством представленного программного комплекса является его полная формализация, позволяющая 
имитировать надежность функционирования любых технологических систем при имеющейся структурной схеме и 
параметрах элементов, входящих в эту схему. 

Проверка предполагаемого подхода по оценке риска, связанного с надежностью исследуемого объекта, про- 
водили на основе имитационного моделирования надежности функционирования универсального измельчителя ком- 
понентов комбикормов типа У17-УКИ. 

Общий вид измельчителя показан на рис. 2, а его структурная схема имитационного моделирования надежно- 
сти представлена на рис. 3. 





Рис. 3 Структурная схема измельчителя У17-УКИ : 1 — бункер для зерна с заслонкой; 2,3 — подшипник в сборе; 
4,5 — уплотнение подшипника; 6 — шкив измельчающего барабана в сборе (со шпонкой и стопорным винтом); 7 — вал 
ротора; 8 — барабан ротора со шпонкой; 9-12 — пальцы молотков; 13 — набор молотков (24 шт.); 14 — роторная втулка барабана 
с винтом; 15 — нож; 16 — терка; 17 — ремень приводной; 18 — электродвигатель трехфазный,; 19, 20 — конденсатор электроли- 
тический; 21 — кнопочная станция; 22 — пускатель; 23 — кабель питания; 24 — шкив электродвигателя в сборе; 25 — загрузочное 
устройство для овощей и зелени; 26 — корпус в сборе 


В таблице 1 представлены фрагменты параметров законов распределения наработок на отказ То времени вос- 
становления Тв и первоначальной стоимости элементов [7]. 





























Таблица 1 
Пример задания исходных данных для имитационного моделирования 
п/п Наименование объекта То, ч Т„, ч Стоимость, руб 
Ср. знач. Закон Ср. знач. Закон 
5 Нож Н Н 
о М=600 в М=0.12 1000 
с =150 с=0.03 
Н Н 
Т — ы 
- ^- 600 М-600 0,12 МО 1000 
с =150 с=0.03 
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Имитационное моделирование проводилось на отрезке времени 10 000 часов при десятикратном повторении. 
Вероятность появления отказов вычислялась на основе данных, полученных в результате имитационного мо- 
делирования работы рассматриваемой системы по формуле статистической вероятности появления отказов 
[8—10]: 
т;(1) 
п(1) 


где т;( 1) — число отказов 1-ого элемента к моменту времени [, п(Е) === общее число отказов всех элементов за проме- 





В()= (5) 


жуток г. Результаты моделирования представлены в таблицах 2 и 3. 


Таблица 2 
Результаты оценки рисков от ненадежности измельчителя У17-УКИ. 
Данные ранжированы по степени риска 







































































№ п/п Название Но Риск, руб 
руб 
13 Комплект молотков 960 227,37 
15 Нож 1000 223,68 
16 Терка 1000 223,68 
18 Электродвигатель 5000 131,58 
24 Шкив ЭД 800 31,58 
7 Ремень 350 18,42 
7 Вал ротора 1200 15,79 
Загрузочное устройство 600 15,79 
22 Пускатель 480 12,63 
19 Конденсатор 1 480 6,31 
20 Конденсатор 2 480 6,31 
21 Кнопочная станция 400 5,26 
10 Палец молотка 3 200 5,26 
23 Кабель электрический 300 3,95 
12 Палец молотка 2 200 2,63 
4 Уплотнение 2 60 0,79 
14 Втулка распорная барабана 24 0,63 





Остальные элементы из схемы на рис. 3 в таблицу не вошли ввиду малого влияния на надежность измельчи- 
теля. 

При моделировании также рассчитывались показатели: коэффициент готовности (К,), общее время восста- 
новления всех отказов (Т,), суммарная стоимость рисков и количество отказов. Усредненные показатели надежности 
функционирования измельчителя У 17-УКИ приведены в таблице 3, а скриншот окна программы показан на рис. 4. 











Таблица 3 
Сводные результаты статистического моделирования надежности 
функционирования измельчителя У 17-УКИ на отрезке времени 10 000 
Средний 
ФФ Общее время Общая стоимость Начальные затраты Количество 
коэффициент 
й восстановления, час рисков, руб (цена машины) отказов 
готовности 
0,99 13,94 69 808 22 394 69 




















Вестник Донского государственного технического университета 2016, №2(85), 84-89 
















































ля наработка Цена Наработка на 'Отказавший Коэффициент Средний Коэффициент Общее время Общая Начальные ы 
}ала отказа отказ системы элемент готовности готовности восстановления стоимость затраты 
> 1455559734 [1000 566.11260816492 Терка 0.999753441563091 0.99464429582515 11.5852757213953 [64414 [23034 
2 1000 32.8813755691235 Нож 0,398836398597102 
3] 24 _ 121,562382345158 Втулка распорная ба... |0.998908004399798 
4 233277577 | 1000 265,573282351458 Нож |0.999563185970554 г 
5 1000 79.3350251050768, Терка 0,998295954062313 
6 1000 347.799520104757 Терка |0. 99978388255152 
7 1000 191.702937712906 Нож 0:999672334230345 
а| [1600 14.1049992042711__ [Комплектмолотков _ [0982938436098776 й — 
3 400 185.77189304908 Кнопочная станция 0:399082850644581 
Н 480 13,1739754154901 Конденсатор 2. 0,983016218992027 
1 000222475 1000 142.41667073125 Терка 0,399258084370303 
1 | 1000 358.588833853232 Нож 0,999626243384744 
1 1000 69.2749150117348 — |Терка 0,997387713741645 
1 |350 20.943955490875 Ремень, [0.993346198937352 
1] 1600 249,680060965233 — (Комплектмолотков _ [0999291440433789 ] | 
1 1000 137.007969545227 Нож 0,399191859751232 
1 [1000 22.9685960597711 Терка 0.993339490530719 
1 284179072 800 226.062594889352. Шкив ЭД 0.398826294050656 
1 480 217.82662706023 Конденсатор 1 0,999247674076243 
2 [1000 183.761954959067 _ Терка 0,99951701537872 
2 [1000 16,1894311140565 Нож 0,391415660094891 
2 [1600 —189.340812872656 _ Комплектмолотов — |0,999129060924754 
2 338700075 | 1000 288.06217625662 Нож 0.999550539206212 
2 1000 303.548688676104 Терка |0,999624032117565 
2| 800 43.2339085965866 Барабан ротора 0.995691488626397 | 
2 4000 25,3059216370475 Корпус в сборе. 0,382756952091007 * 
‹ === кет еее 6 мае Нее Г р 
| Экспортв Ехсе! ] С Построить гистрограмму отказов ] [ Построить гистограмму восстановлений ] [ График средней наработки на отказ ] 




















Рис. 4. Скриншот работы программы 


Проанализируем полученные результаты. Учитывая, что срок гарантийной работы измельчителя составляет 
10 000 часов, то при работе в течение 7 часов в смену 305 дней в году срок службы машины составит 4,68 года. Исхо- 
дя из этих рассуждений и, судя по показателям в табл. 3, машина довольно надежная. Общая стоимость рисков 69 808 
руб. при стоимости машины 22 394 руб., что вполне согласуется с общепринятыми положениями, что эксплуатацион- 
ные затраты обычно в 3—4 раза больше первоначальных. 

Риски рассчитывались, исходя из предположения, что они определяются только стоимостью заменяемой де- 
тали. На самом деле в риски должны включаться стоимость поиска детали и оплата труда рабочего. В среднем это 
составляет 30% от стоимости детали. Стало быть, общая сумма рисков будет 90 750 руб. Посчитаем также упущенную 
выгоду по зерну. Производительность измельчителя по зерну составляет 200 кг/час. За 13,94 часа (общее время вос- 
становления) должно быть измельчено 0,2х13,94=2,78 т. Стоимость одной тонны зерна в настоящее время составляет 
порядка 7 тысяч рублей. Тогда суммарная упущенная выгода составит 110 210 руб. 

Как показывает опыт эксплуатации таких измельчителей, одним из слабых их мест являются молотки. Из ка- 
талога был выбран комплект молотков с наиболее высокой наработкой на отказ. Стоимость его составляла 1 600 руб. 
(против 960 руб. у прототипа). Далее была проведена оценка рисков по той же программе и методике. Результаты 
представлены в таблице 4. 

Таблица 4 
Сводные результаты моделирования надежности функционирования 
измельчителя У17-УКИ с молотками повышенной надежности 








Средний Общее время 
Общая стоимость | Начальные затраты 
коэффициент восстановления, Количество отказов 
рисков, руб. (цена машины) 
готовности час 
0,995 11,58 64 414 23 034 54 























Тогда общая стоимость рисков улучшенной машины с учетом упущенной выгоды 83 746,2 руб., что на 26 463 
руб. меньше, чем было. Таким образом, видим, что увеличение первоначальной стоимости машины на 640 руб. дает 
повышение эффектности на 26 463 руб., т.е. примерно в 30 раз больше, чем вложено. 


Выводы: 

1) Разработана система оценки рисков от ненадежности машин в виде упущенной выгоды. 

2) Метод основывается на принципах статистического моделирования, является полностью формализован- 
ным и может применяться для любых систем и машин. 

3) Проверка метода на основе имитационного моделирования надежности функционирования измельчителя 
компонентов комбикормов типа У17-УКИ показала хорошую сходимость реальных результатов с данными моделиро- 


вания. 
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